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Abstract 


Particles in rotating saddle potentials exhibit precessional motion 
which, up to now, has been explained by explicit computation. We 
show that this precession is due to a hidden Coriolis-like force which, 
unlike the standard Coriolis force, is present in the inertial frame. We 
do so by finding a hodograph-like “guiding center” transformation 
using the method of normal form. 

1 Introduction and background 

We consider the motion of a particle in the rotating saddle potential in the 
plane. The “spinning” potential whose graph is obtained by rotating the 
graph of a hxed potential Uo{x) = Uo{xi,X2) with angnlar velocity oj is 


U{x,t) = Uo{R ^x), X = {xi,X2) G 


where 
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is the counterclockwise rotation. If Uo is a saddle with equal principal cur¬ 
vatures, then 

Uo{x) = -{xl-xl) 

without the loss of generality, and the equations of motion x = —VU take 
the form 

X -|- S{ujt)x = 0, X G (1) 


where 


^(r) 


cos 2r sin 2r 
sin 2r — cos 2r 


r = ut. 


These equations describe the linearized motion of a particle sliding without 
friction on a rotating saddle surface (with equal and opposite principal cur¬ 
vatures) in the presence of gravity. Figured] It is a surprising fact, known for 
almost a century [H E] IS] S] E] |6] that the equilibrium position of the particle 
becomes stable if the surface rotates around the vertical axis sufficiently fast 
(a heuristic explanation of this effect can be found in [7], and is also given 
below). Numerical experiments show another puzzling effect: the Foucault- 
like prograde precession happening in the inertial frame [T] E] |9], Figure |2j 



Figure 1: A particle on a rotating saddle surface. 


A similarly surprising phenomenon is the stabilization of a ball rolling 
without slipping on a rotating saddle surface (several demonstrations can be 
found on YouTube) P ITU] . Superhcially, the two effects appear to be the 
same; however, the reasons for stability are entirely different. For the rolling 
ball, the gyroscopic effect, which has no counterpart for a point mass in our 
case, plays the key role. In fact, the rolling ball is stable even if the surface 
is horizontal and flat, HU. The rolling ball is an entirely different system: 
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first, it is a nonholonomic system (see [12] for more details), unlike the one 
considered here, and second, it has more degrees of freedom. 

Returning to the particle in the rotating saddle potential, the force held 
— R{ut)x in ([T|) admits the following nice interpretation. Consider the saddle 
force held {x, —y), and make it time-dependent by rotating each vector coun¬ 
terclockwise with angular velocity 2u. Equation ([T|) describes the motion of 
the unit point mass in this force held. 

We note that S{t) is the rehection in the x-axis followed by a rotation 
by 2r counterclockwise. Equivalently, S is the rehection with respect to the 
line which forms the angle 2r with respect to the x-axis. Thus S is an 
anti-orthogonal matrix, with eigenvalues ±1, and satisfying 

= /, (2) 

the identity matrix. 

2 The main result: a “hodograph” transfor¬ 
mation. 

The following result was stated in [T3], but without the derivation. 

Theorem 1 Given a vector function x : M — )■ consider its “guiding 

center”, or the “hodograph” image 

u = x-^—S{t/e){x-eJx), J = J ^ . (3) 

If x{t) is a solution of (Qj), then u{t) satisfies 

il - —Ju + —u = e*f{u,u,e), (4) 

where f is a function linear in u, u and analytic in e, in a fixed neighborhood 
of e = 0. The guiding center therefore behaves, (ignoring the 0{e^)-terms) 
as a point charge of unit mass in the potential V (n) = y in the magnetic 
field of constant magnitude B = e^/4 perpendicular to the u-plane. 

The hodograph transformation (]3|) uncovers, in particular, a hidden Coriolis- 
like force — (e:^/4) Jft; this seems to be the hrst example when such a force 
arises in an inertial reference frame. 
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Figure |21^C) shows some trajectories of ([I]) with the motion of the guiding 
center superimposed on them. Some motions of the guiding center itself 
are illustrated in Figure [2]^B). It should be noted that the “magnetic” effect 
e^Ju/A is of higher order than the quadratic restoring force —this 
explains why for e small the “petals” become more closely spaced, Figure |2l 



Figure 2: A typical trajectory x; its guiding center n; their superposition. 


Before proceeding with proof, we make a few observations that become 
clear in the course of the proof. 

Remark 1 Contrary to what one might expect, any averaging transforma¬ 
tion must involve x. In other words, transforming just the conhguration 
X alone can not get rid of the time-dependence in the original equation in 
principle. Put differently, the class of contact transformations is insufficiently 
wide to carry out the averaging of our system (Qp. This is explained in the 
Section O 

Remark 2 The formal averaging procedure due to Kapitsa mi described 
in Landau-Lifshitz |T5] is incomplete and may give an incorrect result if one 
is not cautious or not lucky, as pointed out in Section |5l 

Remark 3 It is tempting to study the system in the rotating frame (as had 
been done, Bi). since in such a frame the equations become autonomous, 
exchanging the time dependence for the two inertial forces: the Coriolis force 
and the centrifugal force (in fact, the system in such coordinates is linear 
and thus admits explicit solutions, [D S1 El ca El m). The motions we 
are studying here correspond to a hxed region of 0 G in the rotating 
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reference system this region is not fixed as oj becomes large, which makes 
this frame poorly suited for using normal form. As a side remark, a physical 
manifestation of the dependence of the neighborhood on u is the fact that 
solutions such as ones in Figure [2] undergo rapid rotations around the origin 
when viewed in a rotating frame. 

3 Rotating saddle in physical applications. 

Before giving a proof of our statement, we mention that time-periodic non- 
autonomous equations ([T]), as well as their autonomous version in the rotat¬ 
ing frame, arise in numerous applications across many seemingly unrelated 
branches of classical and modern physics 13 El DE]. In celestial mechan¬ 
ics the rotating saddle equations describe linear stability of the triangular 
Lagrange libration points L 4 and L 5 in the restricted circular three-body 
problem [20l [21]. By this reason, the classical work by Gascheau of 1843 
may be considered as the hrst one that established stability conditions for a 
particle on a rotating saddle [niEiisg. However, it was not until Brouwer, 
one of the authors of the hxed point theorem in topology, considered in 1918 
stability of a heavy particle on a rotating slippery surface that the 

rotating saddle trap per se became an object for investigation. 

Indeed, according to Earnshaw’s theorem an electrostatic potential can¬ 
not have stable equilibria, i.e. minima, since such potentials are harmonic 
functions. The theorem does not apply, however, if the potential depends on 
time; in fact, the 1989 Nobel Prize in physics was awarded to W. Paul [23] 
for his invention of the trap for suspending charged particles in an oscillating 
electric held. Paul’s idea was to stabilize the saddle by “vibrating” the elec¬ 
trostatic held, by analogy with the so-called Stephenson-Kapitsa pendulum 
[131IT 3 I 23 [26l 123 128] in which the upside-down equilibrium is stabilized by 
vibration of the pivot. Brouwer DEI explicitly demonstrated that, instead 
of vibration, the saddle can also be stabilized by rotation of the potential (in 
two dimensions). This ehect is used, e.g., in quantum optics, in the theory 
of rotating radio-frequency ion traps [7] . 

In plasma physics equations ([T]) appear in the modeling of a stellatron - 
a high-current betatron with stellarator helds used for accelerating electron 
beams in helical quadrupole magnetic helds [S] [161 EHl [30] . In atomic physics 
the stable triangular Lagrange points were produced in the Schrodinger- 
Lorentz atomic electron problem by applying a circularly polarized microwave 


5 


field rotating in synchrony with an electron wave packet in a Rydberg atom 
[2U] . This has led to a hrst observation of a non-dispersing Bohr wave packet 
localized near the Lagrange point while circling the atomic nucleus indef¬ 
initely [31]. Recently, the rotating saddle equations ([T|) reappeared in the 
study of conhnement of massless Dirac particles, e.g. electrons in graphene 
[52] , Even stability of a rotating flow of a perfectly conducting ideal fluid in 
an azimuthal magnetic held possesses a mechanical analogy with the stabil¬ 
ity of Lagrange triangular equilibria and, consequently, with the gyroscopic 
stabilization on a rotating saddle [33]. Finally, we note that in mechanical en¬ 
gineering equations ([T]) describe stability of a mass mounted on a non-circular 
weightless rotating shaft subject to a constant axial compression force [HI5^. 

4 Derivation and proofs. 

4.1 The guiding center transformation. 

In an attempt to bring ([1]) to a normal form, let us choose a new variable 
Xi G via 

X = xi + —S{t/e)xi. (5) 


Remark 4 Heuristically, this transformation is suggested by the following rea¬ 
soning: at a fixed position x, the force S{t/£)x rotates in a circle with angular 
velocity 2uj. If an otherwise free particle were subject to such a force, it would 
move in a circle if we subtract the drift, with |5x| being the centripetal force, and 
thus given by |5x| = (2a;)^r, where r is the radius of the circle. With oj = Ifs this 
gives r = e^|5x|/4, explaining the choice (0j. 

Theorem 2 The transformation ^ converts the system ^\) into the form 

2 3 4 

xi — eSJxi -f- —xi -I- —Jxi - Sxi = 0{e^). (6) 

4 4 16 

Remark 5 The coefficient matrices SJ and S' in (|6]) have zero average; this 
may suggest that the corresponding terms can be killed by some transforma¬ 
tion and yield the averaged equation 

X2 + —X2 -^JX2 = 0{e^), 
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a wrong result which looks almost exactly like the correct equation (jl]), ex¬ 
cept for the sign in front of the Coriolis term. This shows that such formal 
averaging of the coefficients is illegal. 


Proof of Theorem [2] is a routine calculation which we give for the sake of 
completeness. Substituting the derivatives of the transformation (j5]) 


X = xi + ^Sxi + ^Sxi 


X = Xi-\r + 2S'i;i -|- Sxi) 


into ([ID yields 


/ -|- —S )xi H— e^Sxi + —Sxi -|- S'! xi -f —Sxi )= 0. 


Note that S satisfies 


5'(r) = -2S{t)J, S"{t) = -AS{t). 


Using the notations r = t/e, ' = we have ^ so that 

5 = e-^S' = -2e-^SJ 

and 

5 = e-^S" = -Ae-^S 

according to ([ID. Substitution into (jHD gives 

+ — eSJxi — S + —S'xi^ = 0; 

canceling two terms as indicated above and using S'^ = I we get 

I + —S 1 ail — eSJxi + —Xi = 0. 
Multiplying by (/ + ^S)-^ = I - we obtain 


xi — e{ I 


i‘) 


SJxi + —il-—S)xi = 0{e^). 


(7) 


( 8 ) 


(9) 


( 10 ) 


Using = I again turns flTOl) into ([6D, completing the proof of Theorem [2l 

Before proceeding with the rest of the normal form reduction (section [TD, 
we pause to make a few heuristic remarks in the next two sections. 
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5 A dead end transformation: Kapitsa’s av¬ 
eraging. 

It is now tempting to make further change of variables of the form 

xi = {I + e^H{t/e))x2, /c > 1, (11) 

to eliminate the 0(e)-terms in ([6|). Perhaps surprisingly, no such transfor¬ 
mation will work, as can be checked by explicit substitution. This leads to a 
somewhat unfortunate but unavoidable conclusion: 

Theorem 3 No transformation of the form IHJl) can eliminate 0{e)-terms 
in In other words, the class of contact transformations is not rich enough 
to do the averaging in our system. 

We conclude: 

In order to average out the time-dependent terms in the rotating saddle 
trap eguation, it is necessary to widen the class of transformations m to 
include non-contact ones. 

Remark 6 Kapitsa described formal averaging of the scalar eguation 

X + a{t/e)f{x) = 0, 

see also m- This heuristic procedure applies to m as well, although our 
eguation m has a matrix instead of the scalar coefficient a. It must be pointed 
out that Kapitsa’s’ heuristic procedure gives only 0(6“^)-terms, but not the 
cubic ones; and it is these terms that are responsible for the “precession” in 
Figure d. 


6 Effective potential: a heuristic derivation. 

For a fixed location xq, the force vector 

Fo{t) = F{xo,t) = -S{t/e)xo (12) 

rotates counterclockwise. If the force were independent of x near xq, then a 
non-drifting particle would move in a circle, with Fo(f), the centripetal force. 


related to the radius vector r = x — Xq via F = —{2u)‘^r = —{e‘^/4:){x — Xq), 
so that 

X = xo + —F = Xo + —Sxq. (13) 

Now so far we pretended that F given by (|T^ does not depend on x. But 
firsj) gives a better approximation than xq, suggesting a better approximation 
for the force: 

Fi{t) = -S = -Sxq - (14) 

showing that in this approximation (improving on flT^ ) Fi{t) travels in a 
shifted circle, so that the average 

Flit) = - 

this is a restoring force (corresponding to the last term on the left-hand side 
in the averaged equations (jl])). 

To repeat the above in a more intuitive way, consider the force F = 
S{t/e){xo + r{t)) changing with t, as r{t) travels in a small circle counter¬ 
clockwise with the angular velocity 2/e. We have 

S{t/e){xo + r{t)) = Sxo + Sr{t). 

The key point is that the last term is constant: 

FS{t/e)r{t) = r(0), (15) 

thus providing the bias mentioned earlier. This constancy is due to the fact 
that the t-dependence in S(t/e)r(t) enters in both S and r and the two de¬ 
pendencies happen to cancel. Indeed, let us temporarily £x S', while allowing 
r to change; then Sr travels clockwise since S is orientation-reversing. But if 
we £x r instead and let S depend on t, then Sr would rotate counterclockwise 
with the same angular velocity. The two rotations thus cancel, explaining 
the constancy of Sr. 

7 Higher order reduction 

As noted in Section [5], we must widen the class of transformations if we are to 
eliminate the leading time-dependent terms in (jH]). To that end we rewrite 


9 



(jSD as a system 


I xi = yi 

j = eSJyi - ^xi - ^yi + ^Sxi + 0{e^), 
or, more compactly as a system in 

Zi = (Aq + eAi + A2 + £^^3 + Ai + Zi, 


where 


A2 


Zl 


Xx 

yi 



Ax = 


0 0 \ 
Q SJ 



(16) 

(17) 

(18) 


Theorem 4 The “hodograph image” 

E^ 

u = Xi ——SJxi 

of solutions of ^ satisfies 

3 2 

il — —Jii + —u = OiE^^), 
4 4 


(19) 


( 20 ) 


i.e. behaves as a particle in the potential V{u) = and subject to a 

magnetic force in the constant magnetic field B = e^/ 4 perpendicular to the 
u-plane. 


Before giving the proof, we note that Theorem 0] implies Theorem [H Indeed, 
from (j5]) we have 

xi = X — -^Sx + 0 {e^); 

substituting this into flT^ gives 

u = x — —S{x — eJx) + 0 {e‘^), (21) 

with the conclusion that this transformation converts the original equation 
([1]) into fl20l) . This holds even if the 0(£^)-term in fl2ll) is deleted, since it 
affects only the 0(£^)-terms in the equation for u. This proves Theorem [1] 
modulo Theorem m 
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Remark 7 The magnetic term in is of higher order than the restoring 
term. This is confirmed by the numerical evidence given in Figure [B' for 
small e the precession is slow compared to the oscillations. 

Proof of Theorem 01 This theorem can be verified by a direct calculation, 
simply by substituting fET]) into fl2U]) and using the governing equation x = 
—Sx. This, however, would give no hint on how fl21l) was discovered, and we 
proceed with a normal form reduction which led to flT^ (and thus to fl?I]) b 
We wish to eliminate time-dependence of the coefficients of our system 
flTB]) up to and including order We do so by the standard normal form 
argument. 

Averaging the (9(e:)— term. We seek to kill time-dependence in the 0(e)- 
term in 0161) via the change of variables 

Zi = (I + €‘^Ti)z 2 , Ti = Ti(t/e), (22) 

where Ti (r) is a periodic 4x4 matrix function of period vr. Substituting this 
into (1T6|) and using 

(I + = 1- e^Ti + e^Tl + 0(e^), 


we obtain 

Z 2 — (Bq + eBi + B 2 + e^B^ + s^B^ + 0(e^))z2, 

(23) 

where 

-Bq = Aq, Bi = Ai — T(, 

(24) 

and 

B 2 = A 2 + [Aq, Ti] 

53 = A3 + [Ai,Ti]+TiT( 

(25) 


54 = A 4 + [A 2 , Ti] — Ti [Aq, Ti] 



with brackets denoting commutator of matrices. Note that according to our 
notation T' = e:“^T, so that T' = 0(1). By setting 


^■ = -5(0 s) P6) 

we get 

Bi = Ai — T[ = 0, 
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as follows from (1T5]) and Q- Substituting into we compute 


B2 


1 

4 


0 2 ^ \ 

/ 0 y) ’ 



0 0 \ 

0 J y) ’ 


B, 



summarizing, our equation becomes 


Z 2 — {Bq + £^-82 + £^B^ + + 0 (£^))z 2 - (28) 


Averaging of the term. We now eliminate t from the B 2 term in fl28|) 
by seeking the transformation 

Z 2 = (/ + £^T 2 )z 3 , T 2 = T 2 {t/£), (29) 

where T 2 (r) is a matrix function periodic in r of period 27r. Substitution of 
into fl 2 ^ gives the new system 


is = M2Z3 


where 

M 2 = (/ + + £^T2) - (/ + £^T2)-h^T^. (30) 

Note that we used the fact that T 2 = Multiplying out fl30|) and 

collecting the like powers of £ we obtain 

M 2 = Aq + £^{B 2 — T 2 ) + £^ (-83 + [Aq,T 2 ]) + + 0{e^)] (31) 


note that the e^-term was unaffected by the transformation. To kill the 
f-dependence in the e^Merm we choose T 2 so as to turn B 2 — T' into the 
average of B 2 : 

= (32) 

This condition, along with the requirement of periodicity, dictates the choice 


T 2 


1/0 SJ \ 

4 V 0 0 y) • 


(33) 


Substituting this into fITT]) yields 

M2 = Aq + £^ B2 + B3 + + 0(^£^), 

where we used the fact that B 3 + [Ao,T 2 ] = B 3 , since T 2 commutes with Aq. 
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Reduction of the term. Note that the cubic term turned out to be 
time-independent, and thus we need to average the quartic term. To that 
end we subject the system 

is = M2 2:3 

to the transformation 

Z 3 = (/ + 

with the periodic matrix function T 4 chosen so as to kill time dependence in 
The matrix 

Ms = (/ + e^Ts)-^M2{I + e^Ts) - (/ + 
of the transformed system differs from M 2 only in the terms starting with e^: 

Ms = M2 - e^Ti 


and thus we must choose T 4 so as to kill the time-dependence in the coefficient 
of 


Ti = B^ 



0 0 \ 

^ 0 y) ’ 


or 



0 0 \ 
SJ Q )' 


( 34 ) 


Denoting = w, we obtain the averaged system 

w = (do + B 2 + e^B^ + 0 (£^)) w, 


( 35 ) 


or, explicitly. 


d f u 
dt\v 




( 36 ) 


It follows that u satisfies 


3 2 

u-^Ju + ^u = 0{e^); ( 37 ) 

^We use the subscript 4 for consistency, noting that T 3 = 0, i.e. that the identity 
transformation is needed for the cubic terms. 
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indeed, according to the first equation in 

ii = V + 0{e^)] 

differentiating this by t gives 

u = i) + 0{e'^) 

- note the drop in the power of e due to differentiation (recall that d/dt = 
e~^d/dT). Substituting v from the second equation in fl5^ results in fl57D . 

It remains to hnd the explicit form for the averaging transformation 

zi = {I + e^T^){I + e^T2){I + e^Ti)w. 

Expanding the product in the powers of e, we write the transformation as 

I + e^Ti + e^T2 + 0{e^)] (38) 

substituting the expressions for Ti and T 2 (see fl26|) and (jSS])) and reading off 
the hrst component, we obtain ([19]), as claimed in the statement of Theorem 

m 0 

8 Conclusion 

We showed that the rapid rotation of the symmetric saddle potential creates 
a weak Lorentz-like, or a Coriolis-like force, in addition to an effective sta¬ 
bilizing potential - all in the inertial frame. As a result, the particle in the 
rotating saddle exhibits, in addition to oscillations caused by effective restor¬ 
ing force, a slow prograde precession in the inertial frame caused by this 
pseudo-Coriolis effect. By Ending a hodograph-like “guiding center” trans¬ 
formation using the method of normal form, we found the efiective equations 
of this precession that coincide with the equations of the Foucault’s pendu¬ 
lum [35]. Interpretation of the unconventional Coriolis-like force arising in 
the inertial frame in the spirit of the geometric magnetism [3^ [371138] is an 
open problem and would be welcome. 
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